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We analytically work out the perturbation Ap induced by the Kehagias-Sfetsos (KS) 
space-time solution of the Hoi'ava-Lifshitz (HL) modified gravity at long distances on the 
two-body range p for a pair of test particles A and B orbiting the same mass M . We apply 
our results to the most recently obtained range-residuals &p for some planets of the solar 
system (Mercury, Mars, Saturn) ranged from the Earth to effectively constrain the dimen- 
sionsless KS parameter i/jq for the Sun. Wc obtain i/i® > 7.2 X 10"-'^'' (Mercury), i/iq" > 
9 X 10"^-^ (Mars), i/;® > 1.7 X lO"'^^ (Saturn). Such lower bounds are tighter than other 
ones existing in literature by several orders of magnitude. We also preliminarily obtain 
i/jj > 8 X 10"^" for the system constituted by the S2 star orbiting the Supcrmassive 
Black Hole (SBH) in the center of the Galaxy. 

Keywords: Experimental tests of gravitational theories. Modified theories of gravity. 
Celestial mechanics, Ephemerides, almanacs, and calendars. Orbit determination and 
improvement 
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1. Introduction 

Since their publication, the Hofava's seminal paper on a four-dimensional theory 
of gravity which is power-counting renormalizable and, hence, can be considered a 
candidate for the short-range (UV) completion of the General Theory of Relativity 
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2 L. lorio & M. L. Ruggiero 

(GTR), have raised a lot of interest and stimulated an intense research effort on 
many topics of the theory and its modifications (a brief review can be found in 
the recent paper by SotiriorPl). Hoi^ava's theory admits Lifshitz's scale invariance: 
r — >■ 6r, t bH, and, after this, it is referred to as Hofava-Lifshitz (HL) grav- 
ity. Actually, it has anistropic scaling in the UV domain, since it is (7 = 3, while 
relativistic scaling with q = l \s recovered at large distances (IR). 

Among the other issue s , sphe rically symmetric solutions in HL gravity have 
been investigated in detailJ^EEEl^ ajgo in five dimcnsionJ^ (as for slowly rotating 
solutions, see Rcg.[9| [T0|) : in particular, Kchagias and Sfctsoim obtained an asymp- 
totically flat and static spherically symmetric solution that can be considered the 
analog of Schwarzschild solution in GTR. An open issue in HL theory pertains the 
role of matter and its coupling to gravity, which has not been clarified in full details. 
As a consequence, the motion of free particles in HL gravity is not trivial: since the 
fact that particles move alon g geodesies i s not granted, it is possible to expect de- 
viations from geodesic motior 12 | 13 | 14 |T5l jjowever, many authors (sec e.g. Ref. [16] 



and references therein) assumed that free particles followed geodesies of KS met- 
ric to focus on many issues HL gravity, such as gravitational lensing, quasi-normal 
modes, accretion disks and so on (moreover, an accurate analysis of KS geodesies 
can be found in Ref. [T7|) . In other words, they used KS solution as a toy model to 
study some fundamental aspects of the theory. 

In this paper, starting from the same assumption, we aim at con straining the 
dimcnsionsless parameter -^o of KS metric. In previous paperJ^ ^ I ^^ I we obtained 
bounds both from solar system and extra solar system observati ons of orbital mo- 
tions; other constraints were derived by studying light deflectioJISl and analyzing 
the impact on the classical GTR testa^. Here we show that tighter results can be 
obtained by considering the perturbation Ap induced by the KS solution on the 
two-body range p for a pair of test particles A and B orbiting the same central 
body of mass M . 

The plan of the work is as follows: in Section [51 after a brief introduction to 
the KS solution in the context of HL gravity, we will analytically work out the 
perturbation Ap induced by the KS solution on the two-body range p, which is a 
very accurate, direct and unambiguous observable in solar system planetary stud- 
ies. We will try to make such a part as more self-consistent as possible, in view 
of a readership which may not be fully acquainted with the subtleties of celestial 
mechanics applied to fundamental physics, an endeavor requiring interdisciplinary 
knowledge across different fields. In this way the reader has the possibility of fol- 
lowing autonomously the future developments of such kind of investigations, and to 
apply the present approach to other exotic long-range modified models of gravity as 
well. In Section [3| we compare our theoretical predictions to the range residuals 5p 
constructed with the latest planetary ephemerides applied to recent, accurate data 
sets of some planets of the solar system. Section [4| is devoted to summarizing our 
findings. 
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Analytical computation of the two-body range perturbation 

The starting point to obtain the KS solution is the action of the theory, given bj0 

In eq. (HJ, g is the determinant of the metric tensor g-ys, Kij and C*^ are the second 
fundamental form and the Cotton tensor, respectively, whose expressions are 

= 2^ [k, - V,iV, - V,iV,) , (2) 

and 

C'l = e*«Vfc „ , (3) 

in terms of time derivative of the metric of the three-dimensional spatial slices gij , 
the lapse N and the shift iV*, according to the standard ADM formalism (see e.g. 
Ref.El Ch. 21), while and i?^^) 

are the Ricci tensor and scalar of the three- 
geometry. The parameters k,X,w are coupling constants, while A^ is proportional 
to the cosmological constant, anjH the constant v has been introduced to modify 
the original HL action (see below). It is interesting to point out that standard 
GTR is recovered if the running parameter A has a particular value, i.e. if A = 1 
is an IR fixed pointEl Since the natural IR vacuum of the original HL theory is 
anti-de Sitter, in order to investigate the existence of limits of the theory with a 
Minkowski vacuum, Kehagias and Sfetsos deformed the theory with the term i/'^R^^^ 
and, furthermore, they considered the limit Aw 0. In this limit, and looking for 
a statically symmetric solution in the form 

idsr^N{rncdtr-^j^-r^d^)\ (4) 

with ((i$)^ = (d^?)'^-|-sin^ 9{dip)'^ , from ([J) it is possible to obtain the corresponding 
Lagrangian: 



i^v^ N 
~ 8(1-3A) V7 



(2A - 1)^^^^^ - 2A^/' + - 2^(1 - / - rf) 



(5) 



^Latin indices run from 1 to 3, and refer to space components, while Greek run from to 3, and 
refer to spacetime components. 

''Notice that Kehagias and Sfetsos used the symbol fi instead of v; in the present paper /i denotes 
the gravitational parameter of the central body which is the product of the Newtonian gravitational 
constant by its mass. See below, after eq. (0. 

■^However, it is important to stress that even though for A = 1 the action reduces to the one of 
GTR, the theory possesses an extra degree of freedom which becomes strongly coupled when A 
goes to 1 (see e.g. Ref. '22' and references therein). 



April 18, 2011 0:50 WSPC/INSTRUCTION 

HL'range'ML"rev'02032011 



FILE 



4 L. lorio & M. L. Ruggiero 

where ^/^ — 8i'^(3A — 1)/k^. As for the A = 1 case, we get that ijj — IGi/^/k^. By 
solving the equations of motion, the asymptotically flat KS solution turns out to be 



{dsf = e"('')(c(ii)' - e^^'^Hdrf - r\d<i>f, (6) 

where 

and fi = GM is the product of the Newtonian constant of gravitation G times the 
mass M of the central body which acts as source of the gravitational field0 

As we showed in Ref. [THl in typical Solar System conditions, the metric ([7]) 
can be written in the form "Schwarzschild plus corrections". To this end, we set 
f/; = "° , where ipo is a dimensionless parameter; then, we can see that in the limit 



V^o»4(^)'«4xl0-4 (8) 
for M = Mq and r « 1 AU, we obtain 



e . 2 

c r 



As a consequence, on introducing the isotropic radial coordinate r, the metric 
^ becomes 

(10) 

In the following, for the sake of simplicity, we will re-name r as r. 

We point out that, even if here we are studying the KS solution obtained in the 
context of HL gravity, our results can be applied as well to an arbitrary perturbation 
of the Schwarzschild solution in the form e"^''^ « 1 — ^ + ^ where e is a suitable 
perturbation parameter. Also, for the sake of clarity, we stress that in performing 
our calculations we do assume that particles move along geodesic, which is not, in 
general, granted in HL gravity, where also Lorentz invariance does not always hold. 
In summary, we use a relativistic post-Newtonian approach on the KS metric and 
the constraints that we obtain consistently apply to the metric parameters, rather 
than to the HL theory which, as we have seen, shows some inconsistencies. 

From (jiop . in addition to the usual Newtonian and general relativistic 
Schwarzschildian terms, it is possible to obtain the corresponding KS acceleration. 



•^Notice that Kehagias and Sfetsos used the symbol uj instead of we use the latter since the 
symbol ui might be confused with the pericenter of the orbit of a test particle, which is one of the 
standard Kcplerian orbital elements: see the discussion after eq. I|15|l . 
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which reads (see also Ref. [25)) 



A-tflo 



(11) 



where R is the unit vector radially directed towards the test particle [see below, 
eq. (Uni)], while vn is the component of its velocity v along the radial direction (see 
below, eq. HI])). Notice that cq. ([TT|) is only valid at post-Newtonian order: it can be 
viewed as a small correction to the main Newtonian monopole term, to be treated 
with, e.g., the standard perturbative method by Gauss [see below, eq. p^ ]. To this 
aim, let us notice that it usually refers to the Newtonian-Keplerian orbit assumed 
as unperturbed, reference trajectory. Ii i princ iple, it is possible to assume as ref- 
erence path a fully post-Newtonian onJ^HlSI^ and work out the effects of a given 
small extra-acceleration like eq. (|lll) with respect to it according to the perturbative 
scheme set up by the authors of Refs. [24l [25l which is a general relativistic general- 
ization of another standard perturbative approach based on the planetary Lagrange 
cquation^^. Actually, it is, in practice, useless since the only addition with respect 
to the orbital effects like, e.g., the precession of the pericenter, resulting from the 
standard scenario would consist of further, small mixed GTR-perturbation orbital 
effects. To be more explicit, it can be shown that further acceleration terms of the 
order o£| fj,^c~^r~^ipo^ ^^c~^''r~^^(^^ would appear. The bounds on Vo that 
we will obtain in Section |3] will show, a posteriori, that they are, in fact, completely 
negligible in such a way that their inclusion would only make the calculation more 
cumbersome, without substantially affecting the constraints on -00 ■ 

In order to perturbativcly work out the orbital effects of eq. (fTT|) . it must, first, 
be evaluated onto the unperturbed Keplerian ellipsi^H 



1 + e cos / 



where a and e are the semimajor axis and the eccentricity, respectively, character- 
izing its size and shapcj^. while / is the true anomaly reckoning the instantaneous 
position of the test particle with respect to the point of closest approach to the 
central body, generally dubbed pericenter. To this aim, cq. (fTTj) has to be projected 
onto the radial (i?), transverse (T) and normal (TV) directions of an orthonormal 
frame co-moving with the orbitep^. It is spanned by the unit vectors {R, T, N} 
whose components with respect to a locally inertial frame {i, j, k} centered in M 



•^They come from the term — (c^/2)/iij./ioo /si « = 1,2,3 in the post-Newtonian equations of 
motioiI^H. 

*It is < e < 1; circular orbits correspond to e = 0. 
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6 L. lorio & M. L. Ruggiero 



' cos cos It — cos/sinflsinw, 
R — ( sinficosM + cos/cosfisinw, 



(13) 



, sin/ sin u, 



— cos Q, sin u — cos / sin Q, cos 



u, 



T 



sin D, sin u + cos / cos CI cos u, 



(14) 



, sin/ cos u, 



N 



' sin /sin 57, 



sin/ cos r2. 



(15) 



, cos /. 



In eq. p^ -cg. (|15l) / and arc the inclination of the orbital plane with respect to 
the reference {x, y} plane and the longitude of the ascending node, respectively: 
is an angle in the reference {x, y} plane delimited by the reference x axis, usually 
taken coincident with the mean equinox at the epoch J2000.0 in planetary data 
analyses, and the line of the nodes, i.e. the intersection between the reference {a;, y} 
plane and the orbital plane. The angle u = cj + / is the argument of latitude in 
which the angle uj is the argument of pericenter: uj is an angle in the orbital plane 
reckoning the location of the point of closest approach to A/ with respect to the line 
of the nodes. Basically, /, 51, uj can be regarded as the three Euler angles determi ning 
the orientation of a rigid body, i.e. the unperturbed orbil[f|, in the inertial space^S 
By recalling that the R — T — N components of the unperturbed Keplerian test 
particle's velocity ar 

f nae sin f 

< VT = (16) 
. wjv = 0, 



sin the Keplerian two-body problem the ellipse neither changes its shape nor its size, i.e., the 
semimajor axis a and the eccentricity e remain constant. Moreover, it keeps its spatial orientation 
unchanged, i.e., I,Q,uj do not change either. 



April 18, 2011 0:50 WSPC/INSTRUCTION 

HL'rangc'ML'rcv'02032011 



FILE 



Phenomenological Constraints on the Kehagias-Sfetsos Solution 7 



where n = ^ [J^l o? is the unperturbed Keplerian mean motion, it turns out that the 
R~ T — N components of eq. (fTTj) are 

Ar. = iS^S^ {4 (1 - e^) + ^ [(1 - 4e2) + e (2 eos / + he cos 2/)] } , 
At = 



lOe^i (l+ecos/)°sin/ 
i/'oc»ae(l-e^)B ' 



An=0. 

(17) 

Notice that, in the hmit e — > 0, it turns out that, contrary to Aj^, At — > 0: moreover, 
vji — > 0, while vt does not vanish. Actually, there is nothing strange or pathological 
in that: just as in Newtonian case, the acceleration experienced by the test 
particle is entirely radial, and its velocity vector is tangential to the circular tra- 
jectory, changing only its direction. An inspection of the content within the curly 
brackets of Ar in eq. (|T7| shows that the second term is, usually, several orders of 
magnitude smaller that the first one, which is of the order of unity: indeed, in the 
case of the Sun and, say. Mercury it is 

Mo 



= 2 X 10" 



(18) 



Moreover, in the limit e — > 0, valid for circular orbits, as we already noticed. At 
vanishes, while Afi reduces to a constant term depending on the test particle's 
orbital radius as to the leading order. 

At this point, the perturbations of all the six Keplerian orbital elementi^, 
including also the mean anomaly Ai , have to be computed by means of the standard 
Gauss perturbative equations^ 



/ da 

dt 



de 
dt 



dJ 
dt 



[eARsinf + ATi^)] , 



VT 



-^{ARsinf + AT [cos/ + i(l-^)]} 



d.n 

dt 



dip 
dt 



-^^An {^) cosu, 



Am (-] sinti, 



- Vk 



-ArcosJ + At 1 



sin/ 



-cos/f , 



dM 
dt 



(19) 



In eq. (fT9|) p = a(l — e^) is the semi-latus rectum. It turns out that by inserting eq. 



''It also holds in GTR since, as it is well known, the Schwarzschild acceleration becomes radial for 
a circular orbit. 
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([T7)) in eq. ([T5]l and integrating it by means 

dt=l(-)'^=l=d/ (20) 

from an initial value /g to a generic / allows to obtain exact expressions for the 
shifts Aa, Ae, A/, Ail, Aw, AAi. Since they are extremely cumbersome, we will not 
explicitly display them. The inclination / and the node are not affected since, 
according to eq. ([T7]), = 0. 

The further step consists of putting such formulas in the following expressions 
for the R — T — N perturbations of the position vector r taken from Ref . UHl 

Ai? = (^) Aa-acos/Ae+ ^^|AM, 



< AT = a sin / 



Ae + r(cos/ASl + Aw)+ 2_K/r~~^AAl, (21) 



a(l-e^) 

A = 7'(sin u A/ — cos u sin / Afi) . 

In the unperturbed case, there is only the radial component, i.e. r = rR, with r 
and R given by eq. (|T2)) and eq. ((T3)) . respectively. In the case of eq. ((TT)) . only the 
radial and transverse shifts occur because A7V is induced by A/ and Ail which, 
as already noted, vanish. Also in this case, we do not display the resulting exact 
formulas because of their excessive length and complexity. 

Finally, we are ready to analytically work out the perturbation Ap of the two- 
body range p between two test particles A and B orbiting the same central body. 
Indeed, the range is a direct and accurate observable which is widely used in as- 
tronomical studies either by directly ranging from A = Earth to the surface of the 
target body B or to a spacecraft orbiting it. From ° 

= {rA - rs) ■ (rA - re) , 

(22) 

P p 

to be evaluated onto the unperturbed Keplerian ellipses of the two test particles A 
and B, it follows that, for a generic perturbation, the range shift Ap ia^ 

Ap= (Ar'A-AfB)-/5, (23) 

where 

Afj = ARj Rj + AT,- fj + AN.J Nj, j = A, B. (24) 

Note that, in general, all the components of the perturbations of the position vectors 
of both A and B are needed since, in general, Ri ■ Tj and Ri ■ Nj do not vanish for 
i,j = A, B, i ^ j. In other words, Ap is not simply the difference between Ai?A 
and Ai?B, as it might intuitively be guessed. 



April 18, 2011 0:50 WSPC/INSTRUCTION 

HL'rangc'ML'rcv'02032011 



FILE 



Phenomenological Constraints on the Kehagias-Sfetsos Solution 9 



In order to conveniently plot eq. (|23l) as a function of time, we will us e th e 
following useful approximate equation for the true anomalies of both A and 

J ~M + 2es\nM + fe^ sin2A^ + 

(25) 

+ e3 (i|sin37W - \ sinM) + (^sin47W - iisin27W) ; 
indeed, Ai •n{t — tp), where tp is the time of the passage at pericenter. 



3. Confrontation with the observations 

Our analytical calculations, suitably compared with latest observational resultd 
allow to put tighter constraints on the KS parameter with respect to those 
inferred in Ref. [TSI from the correctiontQ Arij to the standard secular perihelion 
precessions. 

By preliminarily working in the limit e — >■ 0, a naive order-of-magnitude eval- 
uation can be obtained by interpreting Aji in eq. ([TT]) as due to an "effective" 
gravitational parameter 

= M + — ^ (26) 

in the Newtonian monopole. Since A/i = ]l[a) — fi cannot be larger than the accuracy 
with which the gravitational parameter of the central body of the system under 
consideration is known, it turns out 

^0 > (27) 

Note that eq. ^7} tells us that the tightest constraints come from test particles 
orbiting very close to highly massive central bodies. A step further consists of taking 
the average of over an orbital period: in doing that we may neglect the terms 
of order 0{c~^) in eq. ()17p . but we do not set e = 0. The result, exact to 0(c~^), is 

-/ X 2(2 + 3e2)/^4 

^("'^)='^+^0C««3(l„, 2)7/2- (28) 

Thus, 

, ^ 2(2 + 36^)/ 

according to which a high eccentricity contributes to further strengthen the con- 
straints on ipQ. 

By applying eq. ()29p to the eight major planets of the solar system we have the 
bounds summarized in Table [T] They are tighter than those in Table 3-Table 4 of 



'Here uj = Q + uj denotes the longitude of pericenter: it is a "dogleg" angle since, in general, it 
does not lie in a plane.. 
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Ref . [TSl by 1 — 6 orders of magnitude, especially for the outer planets. As expected, 
the tightest constrain comes from Mercury, with 

t/-® > 1.1 X 10"", (30) 

which is about 1.4 times larger than the bound in Table 3 of Ref. [THl The use of 
the less accurate eq. (P7)) would have yielded the weaker constraint 

> 8.8 X 10"^^ (31) 

However, such a back-to-thc-envelope approach cannot substitute a more rigor- 
ous treatment. Thus, in the following we will compare our predicted results for the 
KS range perturbation Ap to the post-fi t ra nge residuals of some planets producec^ 



with the latest solar system ephemerideJ^. As we will see, even tighter constraints 
will result from such a detailed analysis. 

It is interesting to note that if we apply eq. ([29]) to the systenlSH constituted 
by the S2 star (a — 1031.69 au, e = 0.8831) orbiting the Supermassive Black Hole 
(/i, = 4 X lO®/i0, CTp, = 0.5 X lO^/i©) located at the center of the Galaxy we get 

^Ao > 8 X 10-1". (32) 

The bound from eq. p7|) would be 

V'o > 2 X 10-12^ (33) 

i.e. 446 times smaller; this clearly shows how eq. ^7} is inadequate for highly ec- 
centric orbits. In order to correctly interpret such result it should be kept in mind 
that the KS parameter -00 1 contrary to, say, G and c, is not a universal constant of 
Nature; in general, it may vary from a space-time to another being, thus, connected 
with the body which acts as source of the gravitational field. 

Moving to the planetary range residuals, the authors of Ref. [30] processed some 
radiometric data from the Messenger spacecraft during its recent flybys of Mercury 
in 2008-2009 with the INPOPlOa ephemerides; according to Table 1 of Ref. [SQl the 
resulting mean and standard deviation of the produced range residuals are 

(5/9= -0.6 ±1.9 m. (34) 

In order to constrain -0® from such data we will use the analytical results of Section 
[5] to calculate a time series for Ap with the same numerical characteristics. The 
result is depicted in Figure [T] it has been obtained for 

i^f = 7.2 X 10^1", (35) 

and has just 

Ap = 0.1 ± 1.9 m. (36) 



Jit is just the case to recall that GTR was fully modelled, so that such residuals account, in 
principle, for any other un-modelled dynamical effects like the KS perturbations. 
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Mercury range: HL-KS perturbation (analyticalcalculation) 




8.2 



S.4 



8.6 



9 



Fig. 1. Analytically computed Earth-Mercury range perturbation Ap caused by KS for ipQ = 
7.2 X 10"^'' over At = 1 yr. The initial conditions, corresponding to J2000.0 (mean ecliptic and 
equinox at such an epoch), have been retrieved from the NASA-JPL WEB interface HORIZONS. 

We checked our analytical result by simultaneously integrating the equations of 
motion of the Earth and Mercury with and without eq. (jlip for the same initial 
conditions used for producing Figure [T| it turns out that the resulting numerically 
computed KS range signal shows an excellent agreement with the analytical one, 
both qualitatively and quantitatively. Thus, we can reasonably be confident in our 
analytical approach. A value for tp^ smaller than that in eq. would generate a 
signal with a larger standard deviation in eq. (|36p . Thus, we take the result of eq. 
psp as our 1 — a lower bound on i/jq . It is about two orders of magnitude tighter 
than the constraint in Table 3 of Ref. UHl It is also 65 times larger than the naive 
result of Table [T] However, we remark that HL gravity was not modelled at all 
in the INPOPlOa ephemcrides, which, as usual, only took into account the main 
well established, standard Newtonian and Einsteinian effects. Thus, if HL really 
existed in Nature, its signature might have been partially removed by the signal 
in the estimation of the state vectors during the data reduction process. In other 
words, there is, in principle, the possibility that ipQ may actually be smaller than 
eq. (|35p . As for any other putative exotic models of gravity, it would be necessary 
to explicitly include HL in the force models used to process the data and repeat the 
entire fitting procedure with such a modified dynamical theory: however, such a task 
is definitely not a trivial one and would require a dedicated analysis by professional 
teams of astronomers. It is outside the scopes of the present paper. On the contrary, 
the approach outlined here is relatively easy to implement, and can be extended to 
other modified models of gravity as well. 

A similar analysis can be done for Mars. Indeed, its range residuals, built with the 
ephemerides INPOPlOa and the data of the Mars Global Surveyor (MGS) spacecraft 



Figure [5] displays the analytically computed KS range signal over the same time 
interval for 




(5/9 = 0.5 ±1.9 m. 



(37) 



^0^ 



,0 - 



= 9 X 10 



,-12 



(38) 
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It turns out that 



Mars range: HL-KS perturbation (analytical calculatioi^ 




2 4 6 8 

t(yr) 



Fig. 2. Analytically computed Earth-Mars range perturbation Ap caused by KS for 'ip^ = 9 X 
10~^^ over Af = 10 yr. The initial conditions, corresponding to J2000.0 (mean ecliptic and equinox 
at such an epoch), have been retrieved from the NASA-JPL WEB interface HORIZONS. 



Ap=:0.3±2m. (39) 

A numerical integration confirms such features, both qualitatively and quantita- 
tively. The constraint of eq. psp is 64 times larger than the approximate one in 
Table [TJ With respect to Table 3 of Ref. [30l eq. (|38)) is more stringent by three 
orders of magnitude. 

Another planet whose orbit's knowledge has greatly been improved in recent 
years thanks to the continuous radiometric ranging to the Cassini spacecraft is 
Saturn. The range residuals over 3 years (2004-2007) produced with the INPOPlOa 
ephemeridefPni yield 

5p = 0.0 ± 17 m. (40) 
Figure [3] depicts the analytically computed KS range perturbation for 

V'lf = 1.7 X 10^12 : (41) 

it has 

Saturn range; HL-KS perturbation (analyticalcalculation) 

30 
20 
10 

f 

-10 
-20 

4 4.5 5 5.5 6 6.5 7 

Kyrl 

Fig. 3. Analytically computed Earth-Saturn range perturbation Ap caused by KS for -i/;® = 
1.7 X 10^^^ over At = 3 yr. The initial conditions, corresponding to J2000.0 (mean ecliptic and 
equinox at such an epoch), have been retrieved from the NASA-JPL WEB interface HORIZONS. 
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Ap = 0.7 ± 17 m. 



(42) 



Also in this case, a numerical integration fully confirms the characteristics of the 
analytical signal of Figure [31 The constraint of eq. (|^T|) is about four orders of 
magnitude tighter than that in Table [1] and six orders of magnitude better than 
Table 4 in Ref. UHl cfr. with the much smaller deviations in the case of Mercury 
and Mars amounting to about one order of magnitude. It is likely so because the 
radiometric data from the inner planets played a major role in the estimation of 
the Sun's gravitational parameter^. As anticipated in Section [21 the lower bounds 
derived for ijjQ make negligible the mixed terms that we did not included in our 
calculation. 

Finally, we note that the solar system constraints on obtained in the present 
paper are more stringent than those derived by the authors of Ref. [TBI through the 
light deflection, and by the authors of Ref. [501 from the orbital motion of Mercury, 
the deflection of light by the Sun, and the radar echo delays. 



4. Summary and conclusions 

In view of the fact that ranging from the Earth to some major bodies of the solar 
system is one of the most direct and accurate way to determine their orbits, in this 
paper we analytically worked out the effects that the Kehagias-Sfetsos solution of 
the Hofava-Lifshitz modified gravity at long distances have on the two-body range 
p. We successfully checked our analytical results by simultaneously integrating the 
equations of motion of the bodies A and B considered by including such exotic 
dynamical effects as well. 

By comparing our predictions with the range-residuals 6p of Mercury, Mars and 
Saturn produced with the recent ephemeridcs INPOPlOa we have been able to effec- 
tively constrain the dimensionless free parameter i/jq entering the Kehagias-Sfetsos 
equations in the case of the Sun. We obtained f/'tf > 7.2 x 10~^° (Mercury), Vcf ^ 
9 X 10~^^ (Mars), -0^ > 1.7 x 10~^^ (Saturn). Such constraints arc orders of mag- 
nitude better than those existing in literature. 

In principle, one should re-process the entire planetary data set by explic- 
itly modeling the Kehagias-Sfetsos dynamical effects in addition to the standard 
Ncwtonian-Einsteinian ones and solving for a dedicated parameter in order to avoid 
the possibility that they, if unmodeled, may be removed from the signal in the es- 
timation of the initial conditions. However, such a task is very time-consuming and 
can effectively be implemented only by skilful and expert teams of specialists in 
astronomical data processing. On the other hand, the relatively simpler approach 
outlined here can be easily and quickly reproduced, and extended to other exotic 
long-range modified models of gravity as well in order to yield reasonable evaluations 
of the magnitude of the effects of interest. 
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Table 1. 
Lower bounds of the KS 
parameter ipQ from eq. 
S29l with o-^ = 1 X 10^ 
m'^ s~-^, from Table 3 of 
Ref.l30l 
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1.1 X 10" 
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Venus 


1.3 X 10" 
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Earth 


5.1 X 10" 
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Mars 


1.4 X 10" 
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Jupiter 


3.6 X 10" 
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Saturn 


5.8 X 10" 
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Uranus 
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Neptune 
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